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In the space environment, instruments onboard of spacecrafts can be affected by 
displacement damage due to radiation. The differential scattering cross section 
for screened nucleus-nucleus interactions - i.e., including the effects due to 
screened Coulomb nuclear fields -, nuclear stopping powers and non-ionization 
energy losses are treated from about 50keV/nucleon up to relativistic energies. 



1. Introduction 

In the space environment near Earth, low energy particles are, for instance, 
found trapped within the radiation belts and are partially coming from 
the Sun. On the other hand, the energies of Galactic Cosmic Rays (GCR) 
extend up to relativistic range. Protons are the most abundant, but alpha 
particles and heavier nuclei are also present (e.g., see Sections 4.1.2-4.1.2.5 
of Ref. [1]). Abundances and energy spectra of GCRs depend on the po- 
sition inside the solar cavity and are affected by the solar activity. Above 
(30-50) MeV/nucleon, the dominant radiation consists of GCRs. At lower 
energies, from 1 MeV/nucleon up to about 30 MeV/nucleon, one also finds 
the so-called Anomalous Cosmic Rays (ACRs). GCRs can reach Earth's 
magnetosphere and interact with upper layers of the atmosphere. These 
interactions produce secondary particles, like for example protons with (10- 
100) MeV energies, which may - in turn - become trapped within the ra- 
diation belts. In addition, during transient phenomena like solar flares and 
coronal mass ejections, Solar Energetic Particles (SEP) are produced in the 
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energy range from few keV's to GeV's. 

All these energetic particles can inflict permanent damages to onboard 
electronic devices employed in space missions. While passing through mat- 
ter, they can lose energy by Coulomb interactions with electrons (electronic 
energy-loss) and nuclei (nuclear energy-loss) of the material. In particular, 
the nuclear energy-loss - due to screened Coulomb scattering on nuclei of 
the medium - is relevant for the creation of permanent defects inside the 
lattice of the material; thus, for instance, it is mostly responsible for the 
displacement damage which is a cause of degradation of silicon devices. 

The developed model - presented in this article - for screened Coulomb 
elastic scattering up to relativistic energies is included into Geant4 distri- 
bution [2] and is available with Geant4 version 9.4 (December 2010). 



2. Nucleus Nucleus Interactions and Screened Coulomb 
Potentials 

At small distances from the nucleus, the potential energy is a Coulomb 
potential, while - at distances larger than the Bohr radius - the nuclear 
field is screened by the fields of atomic electrons. The interaction between 
two nuclei is usually described in terms of an interatomic Coulomb potential 
(e.g., see Section 2.1.4.1 of Rcf. [1] and Section 4.1 of Rcf. [3]), which is a 
function of the radial distance r between the two nuclei 

zZe 2 

V(r) = — * T (r r ), (1) 

where ez (projectile) and eZ (target) are the charges of the bare nuclei and 
is the interatomic screening function. This latter function depends on 
the reduced radius r r given by 

r r = f, (2) 
ai 

where ai is the so-called screening length (also termed screening radius). In 
the framework of the Thomas-Fermi model of the atom (e.g., see Chapters 1 
and 2 of Ref. [4]) - thus, following the approach of ICRU Report 49 (1993) 
-, a commonly used screening length for z = 1 incoming particles is that 
from Thomas-Fermi (e.g., see Refs. [5, 6]) 



Ctf &o 



(3) 



January 11, 2011 2:42 WSPC - Proceedings Trim Size: 9in x 6in Consolandi'Rancoita 



and - for incoming particles with z > 2 - that introduced by Ziegler, Biersack 
and Littmark (1985) (and termed universal screening length^) 

_ C'tf ap , . 

au — 2,0.23 _|_ _£0.23 ' V 4 J 

where 

h 2 

a = — o 



is the Bohr radius, m is the electron rest mass and 

Ctf = \ (^j ' 0.88534 

is a constant introduced in the Thomas-Fermi model. 

The simple scattering model due to Wentzel [10] - with a single ex- 
ponential screening-function \&i(r r ) {e.g., see Ref. [10] and Equation (21) 
in Ref. [11]} - was repeatedly employed in treating single and multiple 
Coulomb-scattering with screened potentials (e.g., see Ref. [11] - and refe- 
rences therein - for a survey of such a topic and also Refs. [12-15]). The 
resulting elastic differential cross section differs from the Rutherford diffe- 
rential cross section by an additional term - the so-called screening parame- 
ter - which prevents the divergence of the cross section when the angle 9 of 
scattered particles approaches 0°. The screening parameter A Sy M [e.g., see 
Equation (21) of Bcthc (1953)] - as derived by Moliere (1947, 1948) for the 
single Coulomb scattering using a Thomas-Fermi potential - is expressed 13 
as 

2" 



(5) 



where a\ is the screening length - from Eqs. (3, 4) for particles with z = 1 
and z > 2, respectively; a is the fine-structure constant; p (/3c) is the 
momentum (velocity) of the incoming particle undergoing the scattering 
onto a target supposed to be initially at rest; c and H are the speed of light 
and the reduced Planck constant, respectively. When the (relativistic) mass 



a Another screening length commonly used is that from Lindhard and Sharff (1961) 
(e.g., see Ref. [8] ; see also Ref. [9] and references therein): 

_ Ctf an 

^ " (W3 +Z 2/3)l/2- 



b It has to be remarked that the screening radius originally used in Refs. [12, 13] was 
that from Eq. (3). 
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- with corresponding rest mass m - of the incoming particle is much lower 
than the rest mass (M) of the target nucleus, the differential cross section 

- obtained from the Wentzel-Moliere treatment of the single scattering - is: 

da WM (6) fzZe 2 \ 2 1 

(6) 



dfl \p/3c J {2A s .m + 1 - cosf?) 2 

(7) 



zZe 2 \ 2 1 



{*P0cJ [A sM +sm 2 (6/2)Y 

(e.g., see Section 2.3 in Ref. [11] and references therein). Equation (7) dif- 
fers from Rutherford's formula - as already mentioned - for the additional 
term A s ,m to sin 2 (6>/2). The corresponding total cross section {e.g., see 
Equation (25) in Ref. [11]} per nucleus is 



zZe 2x 2 



pfic) A s>m (1 + ^ s ,m) 



(8) 



Thus, for (3 ~ 1 (i.e., at very large p) and with A s _m <C 1, from Eqs. (5, 8) 
one finds that the cross section approaches a constant: 



ff WM 



2zZe 2 a^ 2 



fic J 1.13 + 3.76 x (azZ) 



2 • 



(9) 



In case of a scattering under the action of a central potential (for in- 
stance that due to a screened Coulomb field), when the rest mass of the 
target particle is no longer much larger than the relativistic mass of the inco- 
ming particle, the expression of the differential cross section must properly 
be re-written - in the center of mass system - in terms of an "effective par- 
ticle" with momentum (p'A equal to that of the incoming particle (p' in ) and 
rest mass equal to the relativistic reduced mass 

_ mM 
Mrcl " Mi? 

where is the invariant mass; m and M are the rest masses of the inco- 
ming and target particles, respectively (e.g., see Refs. [15-17] and references 
therein). The "effective particle" velocity is given by: 



1 + 
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Thus, the differential cross section f per unit solid angle of the incoming 
particle results to be given by 



da WM (6') ( zZe 



<KV \2p' in p T cJ [A s + sm 2 (6'/2)] 



with 



-4 S 



2 Pin °I 



1.13 + 3.76 x 



azZ 



(10) 



(11) 



and 0' the scattering angle in the center of mass system. 

The energy T transferred to the recoil target is related to the scattering 
angle as T = T max sin (9'/2) - where T max is the maximum energy which 
can be transferred in the scattering (e.g., see Section 1.5 of Ref. [1]) -, thus, 
assuming an isotropic azimuthal distribution one can re-write Eq. (10) in 
terms of the kinetic energy transferred from the projectile - , i.e., [— T], 
where the negative sign indicates that energy is lost by the projectile - to 
the recoil target as 

rf/ M - ( 7 r^4~) 2 o sm(0') d6' r d<j> 

\2PinPrc) [A s + sm 2 (9'/2)] 2 Jo 

Pin&c) [T max A s + T] 2 d[ ~ n (12) 

Finally, from Eq. (12), the differential cross section with respect to the 
kinetic recoil energy (T) of the target is given by: 

da WM (T) zZe 2 \ 2 T max 



dT \PinPrCj [T max A s + T] 



f By inspection of Eqs. (5, 7, 10, 11), one finds that for (3 r «l the cross section is given 
by Eq. (9). 

c One can show - e.g., see Section 1.5 of Ref. [1] - that the four momentum transfer is 
given by 

t = -2MT. 

Since t is invariant, then the kinetic energy transferred is also invariant. Furthermore, 
since T = T max sin 2 (0'/2), then one finds that 

d[-T] = T max d[- sin 2 (e'/2)] = sin (9') d8> 

(e.g., see Section 1.5 of Ref. [1]). 
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Furthermore, since 



pM 



Pin = 



T 



Mi, 2 

2p 2 M 

"W7 



(14) 



with p and E the momentum and total energy of the incoming particle in 
the laboratory, then one finds 

IE 2 



Therefore, Eq. (13) can be re-written as 



dT 



P 2 Mc* [T max A s +Tf 



(15) 



Equation (15) expresses - as already mentioned - the differential cross sec- 
tion as a function of the (kinetic) energy T achieved by the recoil target. 
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Fig. 1. Nuclear stopping power - in McVcm 2 g 1 - calculated using Eq. (17) in sili- 
con is shown as a function of the kinetic energy per nucleon - from 50 keV/nucleon up 
lOOTeV/nuclcon - for protons, o-particlc and n B-, 12 C-, 28 Si-, 56 Fc-, 115 In-, 208 Pb- 
nuclei. 
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Fig. 2. Nuclear stopping power - in MeVcm 2 g 1 - calculated using Eq. (17) in lead 
is shown as a function of the kinetic energy per nuclcon - from 50 keV/nucleon up 
lOOTeV/nucleon - for protons, a-particles and n B-, 12 C-, 28 Si-, 56 Fe-, 115 In-, 208 Pb- 
nuclei. 

3. Nuclear Stopping Power 

Using Eq. (15) the nuclear stopping power - in MeVcm -1 - is obtained as 



( dx ) nud UA L 



da WM (T) 



dT 



TdT 



= 2 n A ir (zZe 



E 2 



= 2n A Tr(zZe 2 y 



p 2 Mc 4 
E 2 



[A a T„ 



A s + 1 



- 1 + ln 



T] 

'A, 



■dT 



1 



A s 



(16) 



(17) 



p 2 Mc 4 

with ua the number of nuclei (atoms) per unit of volume and, finally, the 
negative sign indicates that the energy is lost by the incoming particle 
(thus, achieved by recoil targets). For energies higher than a few tens of 
keVs, because A s <C 1, Eq. (17) can be re-written as 



(-) 

\ dx J nucl 



E 2 



1 



(18) 



p 2 Mc 4 

It has to be noted that, as the incoming momentum increases to a value 
for which p ~ E, the set of terms - in front of those included in brackets - 
decreases and approaches a constant; while the term \n(l/A s ) increases as 
ln(p) for E > mc 2 ,Mc 2 [e.g., see Eqs. (11, 14)]. Thus, a slight increase of 
the nuclear stopping power with energy is expected because of the decrease 
of the screening parameter with energy. 
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Fig. 3. Universal stopping power (dashed line) as a function of the universal reduced 
energy [Eq. (19)]. The other curves correspond to the dimensionless nuclear stopping 
power obtained from Eq. (17) - for protons, a-particles, 28 Si-, 56 Fe-, 208 Pb-nuclei in 
silicon (Fig. 1) and lead (Fig. 2) absorbers - and divided by the parameter K [Eq. (21)]. 



For instance, in Fig. 1 (Fig. 2) the nuclear stopping power in silicon 
(in lead) - in MeVcm 2 g _1 - is shown as a function of the kinetic energy 
per nucleon - from 50keV/nucleon up 100 TeV/nucleon - for protons, a- 
particles and n B-, 12 C-, 28 Si-, 56 Fe-, 115 In-, 208 Pb-nuclei. 

It has to be remarked that - at very low energies - the Wentzel-Moliere 
nuclear stopping power [Eq. (17)] differs from that obtained by Ziegler, 
Biersack and Littmark (1985) using the so-called universal screening po- 
tential (see also Ref. [18]). However, they have shown (e.g., see Figure 2-18 
in Ref. [7] or, equivalently, Figure 2-18 in Ref. [18]) that different screening 
potentials - including the Bohr potential in which ^\(r r ) is assumed to be 
an exponential function similarly to Wentzel's assumption - result in nu- 
clear stopping powers which exhibit marginal differences for e r .u above 10 
(see also Fig. 3). e rj u is the so-called universal reduced energy expressed as: 

M 

e r,U 



E k , 



(19) 



zZ ( z 0.23 + Z 0.23) \m + Mj 

where is in MeV and the numerical constant is I 
{e.g., see Equation (2-73) of Ref. [7] or Equation (2- 
Section 2.1.4.1 of Ref. [1]}. For instance, in silicon e r u ~ 10 corresponds 
to Ek — 13keV [67 keV/ nucleon] for protons [lead nuclei]. Ziegler, Biersack 
and Littmark (1985) provided a general expression for the nuclear stopping 



: 32.536xl0 3 MeV- 1 
of Ref. [18], see also 
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Fig. 4. Variation (in percentage) of nuclear stopping powers - calculated with Eq. (17) 
for energies from 50 keV/nucleon up to 100 MeV/nucleon - with respect to ICRU tabu- 
lated values [3] as a function of the kinetic energy in MeV/nucleon, for protons (top) and 
ct-particles (bottom) traversing amorphous carbon, silicon, germanium and lead media. 



power (e.g., see Section 2.1.4.1 of Rcf. [1]), i.e., 

— = K<K(e rU ) [MeV/cm], 
dx I „„„, 



(20) 
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Fig. 5. Variation (in percentage) of nuclear stopping powers - calculated with Eq. (17) 
with the expressions (22, 23) for the screening parameter - with respect to ICRU tab- 
ulated values [3] as a function of the kinetic energy in MeV/nucleon, for o-particlcs 
traversing amorphous carbon, silicon, germanium and lead media for energies from 
50kcV/nucleon up to 100 MeV/nucleon. 
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with 



K ~ 5.1053 xlO 3 



p zZ 



[MeV/cm] 



(21) 



A(l + M/m) (z - 23 + Z - 23 ) 



where p and A are the density and atomic weight of the target medium, 
respectively; 9i(e r ^u) - the so-called (universal) reduced nuclear stopping 
power [termed, also, (universal) scaled nuclear stopping power] given in 
Equations (2-89)-(2-90) in Ref. [18] (see also page 80 in Ref. [1]) - is di- 
mensionless. Additionally, the present calculations can be compared with 
values of stopping powers - obtained using the universal screened potential 

- available in SRIM (2008) [19]. Usually an agreement - to better than a few 
percents - is achieved down to about 150keV/nucleon, where - for instance 

- one finds « 5.5 (9.9) % for a-particles (lead ions) in silicon. At large ener- 
gies, the non-relativistic approach due to Zicgler, Biersack and Littmark 
(1985) becomes less appropriate and deviations from stopping powers cal- 
culated by means of the universal screening potential are expected and 
observed for e riU > (1.5-2.5) x 10 4 (e.g., see Fig. 3). 

The non-relativistic approach - based on the universal screening po- 
tential - of Zicgler, Biersack and Littmark (1985) was also used by 
ICRU (1993) to calculate nuclear stopping powers - currently available 
on the web (e.g., see Ref. [20]) - due to protons and a-particles in ma- 
terials. ICRU (1993) used as screening lengths those from Eqs. (3, 4) for 
protons and a-particles, respectively. In Fig. 4, the variation (in percent- 
age) of nuclear stopping powers - calculated with Eq. (17) - with respect 
to ICRU tabulated values [3] is shown as a function of the kinetic energy 
per nucleon (in MeV/nucleon) - for energies from 50keV/nucleon up to 
lOOMeV/nucleon - for protons and a-particles traversing amorphous car- 
bon, silicon, germanium and lead media. The stopping powers for protons 
(a-particles) from Eq. (17) are less than ss 5% larger than those reported 
by ICRU (1993) from 50kcV/nuclcon up to « 8MeV (19 MeV/nucleon) - 
the upper energy corresponds to e rj u wl.6 x 10 4 (2.6 x 10 4 ) for protons in 
an amorphous carbon (a-particles in a silicon) medium -. At larger ener- 
gies the stopping powers from Eq. (17) largely differ from those from ICRU 

- as expected - due to the complete relativistic treatment of the present 
approach. 

The simple screening parameter used so far [Eq. (11)] - derived by 
Molicre (1947) - can be modified by means of a practical correction, i.e., 



( 



h 



) 



2 




azZ 



) 



2 




2 Pin «I 



1.13 + 3.76 x C 



(22) 
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to achieve a better agreement with low energy calculations of Ziegler, Bier- 
sack and Littmark (1985). For instance, for a-particles and heavier ions, 
with 

C = (lOnzZaf 12 (23) 

the stopping powers obtained from Eq. (17) - in which A' s replaces A s - 
differ from the values of SRIM (2008) by less than w 4.7 (3.6) % for a- 
particles (lead ions) in silicon down to about 50 keV/nucleon. With respect 
to the tabulated values of ICRU (1993), the agreement for a-particles is 
usually better than 4% at low energy down to 50kcV/nuclcon (Fig. 5) 
- a 5% agreement is achieved at about 50keV/nuclcon in case of a lead 
medium. At very high energy, the stopping power is slightly affected when 
A' s replaces ^4 S : for example, a) in silicon at lOOTeV/nuclcon the nuclear 
stopping power of a-particles (lead-ions) is decreased by about 0.5 (1.9) % 
and b) in lead it is decreased by about 1.0 (2.7) %. It has to be remarked that 
a more appropriate expression for the screening parameter and a practical 
correction factor may require a further understanding. 

4. Non-Ionizing Energy Loss due to Coulomb Scattering 

A relevant process - which causes permanent damage to the silicon bulk 
structure - is the so-called displacement damage (e.g., see Chapter 4 of 
Ref. [1], Rcfs. [21-23] and references therein). Displacement damage may 
be inflicted when a primary knocked-on atom (PKA) is generated. The in- 
terstitial atom and relative vacancy are termed Frenkel-pair (FP). In turn, 
the displaced atom may have sufficient energy to migrate inside the lat- 
tice and - by further collisions - can displace other atoms as in a colli- 
sion cascade. This displacement process modifies the bulk characteristics 
of the device and causes its degradation. The total number of FPs can be 
estimated calculating the energy density deposited from displacement pro- 
cesses. In turn, this energy density is related to the Non-Ionizing Energy 
Loss (NIEL), i.e., the energy per unit path lost by the incident particle due 
to displacement processes. 

In case of Coulomb scattering on nuclei, the non-ionizing energy-loss can 
be calculated using the Wentzel-Molierc differential cross section [Eq. (15)] 
discussed in Sect. 2, i.e., 

where E is the kinetic energy of the incoming particle, T is the kinetic 
energy transferred to the target atom, L(T) is the fraction of T deposited 
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Fig. 6. Non-ionizing stopping power - in MeVcm 2 g _1 - calculated using Eq. (24) in 
silicon is shown as a function of the kinetic energy per nucleon - from 50 keV / nucleon up 
lOOTeV/nuclcon - for protons, a-particlcs and n B-, 12 C-, 28 Si-, 56 Fe-, 115 In-, 208 Pb- 
nuclei. The threshold energy for displacement is 21 eV in silicon. 
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Fig. 7. Non-ionizing stopping power - in MeV cm 2 g 1 - calculated using Eq. (24) in 
lead is shown as a function of the kinetic energy per nucleon - from 50 keV/nucleon up 
100 TcV/nuclcon - for protons, a-particles and n B-, 12 C-, 28 Si-, 56 Fe-, 115 In-, 208 Pb- 
nuclci. The threshold energy for displacement is 25 eV in lead 



by means of displacement processes. The expression of L{T) - the so-called 
Lindhard partition function - can be found, for instance, in Refs. [24, 25] 
and in Equations (4.94, 4.96) of Section 4.2.1.1 in Ref. [1] (see also referen- 
ces therein). — T L(T) is the so-called damage energy, i.e., the energy 
deposited by a recoil nucleus with kinetic energy T via displacement da- 
mages inside the medium. The integral in Eq. (24) is computed from the 
minimum energy Td - the so-called threshold energy for displacement, i.e., 
that energy necessary to displace the atom from its lattice position - up to 
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Fig. 8. Non-ionizing stopping powers of protons in silicon with energies from 100 kcV 
up to 1 GeV: data for the dashed-dotted and dotted lines are the relativistic and Zicglcr- 
Bicrsack-Littmark screened potential calculations from Jun and collaborators (2003), 
respectively (e.g., see Ref. [26]); data for the dashed line are obtained from Eq. (17) (e.g, 
see Fig. 1). 



the maximum energy T max that can be transferred during a single collision 
process. Td is about 21 eV in silicon (e.g., see Table 1 in Ref. [26] and refe- 
rences therein) and 25 eV in lead (e.g., see Table 22 at page 83 in Ref. [27] 
and references therein). For instance, in Fig. 6 (Fig. 7) the non- ionizing 
energy loss - in MeV cm 2 g -1 - in silicon (lead) is shown as a function of 
the kinetic energy per nucleon - from 50keV/nucleon up 100 TeV/nucleon - 
for protons, a-particles and n B-, 12 C-, 28 Si-, 56 Fe-, 115 In-, 208 Pb-nuclei. As 
already remarked by Boschini and collaborators (2010), at high energy the 
Coulomb NIEL - similarly to the nuclear stopping power - does not de- 
crease with energy as it is found by Ziegler, Biersack and Littmark (1985) 
or in other calculations based on their universal screening potential derived 
in the framework of a non-relativistic treatment of the screened Coulomb 
scattering. 

Furthermore, Jun and collaborators (2003) have already demonstrated 
that a relativistic treatment [24] of Coulomb scattering of protons - with 
kinetic energies above 50 MeV - upon silicon results into a non- ionizing 
energy loss which is larger than that expected from calculations using the 
Zieglcr-Biersack-Littmark screened potential with a universal screening 
length (e.g., see Refs. [7, 25, 26]). The relativistic cross section used for 
treating the Coulomb scattering is the one derived by McKinley and Fes- 
hbach (1948) to describe the scattering of electrons on nuclei (e.g., see 
Section 4.2.1.4 of Ref. [1] and references therein). Seitz and Koehler (1956) 
suggested that - when the mass of the projectile is much lower than the 
target rest-mass (e.g., see Section 13 of Ref. [29] and references therein) 
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- this cross section can also describe - although screening effects are ne- 
glected - the scattering of protons and light nuclei, thus, providing - at 
high energy - a damage cross section which does not decrease with in- 
creasing energy. The data from Jun and collaborators (2003) - for protons 
with energies from 100 keV up to 1 GeV - are shown in Fig. 8: the Ziegler- 
Biersack-Littmark (1985) screened potential was used to treat the Coulomb 
scattering of protons with energies lower than 50 MeV. In the same figure, 
the data obtained using Eq. (17) - e.g., see Fig. 1 - are also shown. There is 
an agreement to better than « 6.5% - achieved at « 1 GeV - between the 
results obtained by Jun and collaborators (2003) and the present calcula- 
tions. 

5. Conclusions 

The treatment of nucleus-nucleus interactions due to relativistic Coulomb 
scatterings with screened potentials - like the present approach based on 
Wentzel-Moliere scattering - allows one to determine both the total and 
differential cross sections, thus, to calculate the resulting nuclear and non- 
ionizing stopping powers. At high energies, the nuclear stopping powers 
exhibit a very slight logarithmic increase with energy. At low energies - i.e., 
above w 50keV/nucleon and up to ~ 8 and 19 MeV/nucleon, for protons 
and a-particles, respectively -, the present results are in agreement to bet- 
ter than w 5% with ICRU tabulated values of stopping powers obtained 
using the universal screening potential for the scattering of protons and 
a-particles in matter. An agreement to better than to (5.5-9.9) % - for in- 
stance, for a silicon medium - is found with SRIM (2008) values down to 
150keV/nucleon. Furthermore, these calculations are also in agreement to 
better than w 6.5% with those obtained by Jun and collaborators (2003) for 
the non-ionizing stopping powers of protons - with energies from 50 MeV 
up to 1 GeV - in silicon. 

Finally, it has to be remarked that - with a simple correction fac- 
tor applied to the screening factor - an agreement to a few percents 
can be achieved with SRIM (2008) stopping powers down to about 
50keV/nucleon. However, a more appropriate expression for the screening 
parameter and a practical correction factor may require a further under- 
standing. 
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